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Abstract 

An ongoing effort to compute the topological susceptibility for the 
f^S SU (3) Yang-Mills theory in the continuum limit with a precison of 

about 2% is reported. The susceptibility is computed by using the def- 
inition of the charge suggested by Neuberger fermions for two values of 
the negative mass parameter s. Finite volume and discretization effects 
are estimated to meet this level of precision. The large statistics re- 
quired has been obtained by using PCs of the INFN-GRID. Simulations 
with larger lattice volumes are necessary in order to better understand- 
ing the continuum limit at small lattice spacing values. 
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1 Introduction 



We continue our numerical study on the lattice of the topological charge dis- 
tribution in the SU(3) Yang-Mills theory adopting the definition suggested 
by Neuberger fermions as discussed in a series of papers pQ-[6]. 

Recent numerical studies on this subject can be found in refs.0 13 [9] 
and ref. [TO]. In the last paper a systematic study at different volumes and 
values of the lattice spacing was performed in order to obtain a reliable 
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determination of the topological susceptibility, the first cumulant of the dis- 
tribution, at the 5% level in the continuum limit at infinite volume. The 
result supports the Witten-Veneziano explanation for the large mass of the 
i- 

The aim of refs [T2] was to look for non-gaussianities in the topologi- 
cal charge distribution of the SU{3) Yang-Mills theory. This is particularly 
challenging since the contribution of the n th cumulant to the charge distri- 
bution is suppressed as V 1 ^ 1 in its asymptotic expansion, even though the 
cumulant itself is a quantity of 0(1) in the infinite volume limit. It is also 
worth noting that in order to search for such very small sub-leading effects it 
is necessary to be sure that all the systematics of the calculation cannot ei- 
ther simulate or hide the effect, and therefore a solid theoretical framework, 
such as the one provided by the topological charge definition suggested from 
Ginsparg- Wilson fermions, is indispensable. In ref. [TU] three different lat- 
tices at the same physical volume of ~ (1.12fm) 4 but with about ten time 
more statistics than in ref. [10], i.e. roughly 3- 10 4 configurations, were stud- 
ied in order to be able to unveil deviations from the Gaussian distribution. 
Significant deviations were found, the second cumulant of the distribution 
gets a value definitively different from zero within errors, and the distribu- 
tion function agrees at that precision level with a first order modification in 
V^ 1 of the Gaussian (Edgeworth expansion). The results clearly disfavour 
the 6 behaviour of the vacuum energy predicted by dilute instanton models, 
while they are compatible with the expectation from the large N c expansion. 

Here we present preliminary results for a measurement on the lattice of the 
topological susceptibility at the 2% level. In order to keep finite volume and 
discretization effects below this level of precision, we have exploited data 
produced in ref. [T2| together with new data generated by additional lat- 
tices as described below. 

All the above challenging Monte Carlo calculation have been made possible 
by important improvements in algorithms which guarantee the reliability 
and the feasibility of high statistics. In particular we have used algorithms 
for zero mode counting with no contamination from quasi zero modes, op- 
timized to run fast on a single processor |13j . The considerable amount of 
computer time needed has been granted to us by the INFN GRID project. It 
allowed us to use the computer resources shared in the scientific Italian net- 
work provided by INFN along these years. We have also taken in advantage 
of the computer resources of the Italian organization COMETA. 



2 



2 Theoretical framework 



In this section we summarize our notation and the necessary theoretical 
framework, for a complete discussion see for example [12J. In the following 
we use the plaquette Wilson action of the SU (3) gauge field. The massless 
lattice Neuberger-Dirac operator D satisfies the Ginsparg- Wilson relation 

mm 

75^ + £>75 = aDj 5 D , (1) 
and the associated topological charge density can be defined as 

a 4 q(x) = -^Ti[ l5 D(x,x)], (2) 

where a = a/(l + s), a is the lattice spacing and s is the negative mass 
parameter. The latter has been fixed in our calculation at the values s = 0.4 
for the study of the volume effects and to s = 0.0 and s = 0.4 for the study 
of the discretization behaviour. The topological charge is obtained from 
the lattice by computing on each gauge configuration the number and the 
chirality of the zero modes of D with the algorithm proposed in Ref. |13j . 
The index v of the Dirac operator 

v = n + — n_ (3) 

is directly related to the topological charge Q 

v = Q = a 4 ^2q(x) . 

X 

In the Euclidean space-time the ground-state energy F{9) is defined as 

e -F(0) = {e i6 Q) (4) 

where, as usual, (...) indicates the path-integral average (our normalization 
is -F(O) = 0). In the large volume regime F(9) is proportional to the size V 
of the system, a direct consequence of the fact that the topological charge 
operator Q is the four-dimensional integral of a local density. The function 
F{9) is related to the probability of finding a gauge field configuration with 
topological charge Q = v by the Fourier transform 

P v = I" M e -"» e -m . (5) 

Large -/V c arguments with N c being the number of colors, suggest that the 
fluctuations of the topological charge are of quantum non-perturbative na- 
ture. The 9 dependence of the vacuum energy is expected at leading order 
in 1/N C , and the normalized cumulants 

1 d 2n 
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which should scale asymptotically as N 2 ~ 2n , have to be determined with a 
non-perturbative computation. The normalized cumulants C n can thus be 
defined as the integrated connected correlation functions of n charge den- 
sities (correlation functions of an odd number of topological charges vanish 
thanks to the invariance of the theory under parity) : 

f .8n 

Cn = y ^ ( q ( Xl )... q (x 2n )r n . (7) 

Xl,..;X2n 

They have an unambiguous finite continuum limit which is independent of 
the details of the regularization [HI [5j |6] . At finite lattice spacing they are 
affected by discretization errors which start at 0(a 2 ). 

The Monte Carlo technique adopted here generates the gauge configu- 
rations with a probability density proportional to e~ SYM , with Sym being 
the chosen discretization of the Yang-Mills action. This algorithm performs 
an importance sampling of the topological charge with the probability dis- 
tribution given in Eq. (JsJ) . A statistical signal for the n th cumulant is then 
obtained when the number of configurations in the sample is high enough to 
be sensitive to terms suppressed as V n ~ l in the asymptotic expansion. For 
instance, the estimators of the first two cumulants 

i=i 

N I N \ 2 

1=1 \ i=l / 

with Vi being the value of the topological charge for a given gauge configura- 
tion and N the total number of configurations, have variances which, up to 
sub-leading corrections, are given by (2<r 4 + a 2 r)/N and (24<r 8 + 72a®r)/N 
respectively being a 2 = VC\ and r = C%IC\. 

3 Lattice data and results 

The properties of the lattices considered and the results obtained for the 
topological susceptibility are reported in Table [TJ 

On the sets Ai, Bo, Bi, B2, B3, Ci and Di the topological charge is com- 
puted with s = 0.4, while it is determined with s = 0.0 on the lattices 
Bo, 0i, £?2 and B3. In order to estimate the magnitude of finite size effects 
we have considered four lattices (Ai, Bi, Ci, Dj)) at the same value of 
the coupling constant corresponding to j3 = 6.0 and at the same value of 
the negative mass parameter s = 0.4. The values of the topological charge 
rescaled with respect to the reference volume of 16 4 are shown in Fig. [Tj It 
is rather clear that for these volumes the results are scattered in a 2% band 
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Lat 


P 


L/a 


r /a 


L[fm] 


N 


(Q' 2 ) 


4 

r oX 


Ai 


6.0 


12 


5.368 


1.12 


34800 


1 627(43") 


0652(421 


Bn 


5.9138 


12 


4.601 


1.30 


10000 


3.271(47) 


0.0707(15) 


Bi 

1 J 1 


6.0 


14 


5.368 


1.30 


30000 


3.097(26) 


0.0669(12) 


Ro 
jJi 


6 0808 


1 6 


6 1 35 


1 30 


1 0000 


? Ql 4.(4.1) 


0630(4 41 


B, 


6.1568 


18 


6.902 


1.30 


10000 


2.843(42) 


0.0615(14) 


Ci 


6.0 


16 


5.368 


1.49 


10000 


5.314(75) 


0.0673(14) 


Di 


6.0 


18 


5.368 


1.60 


10000 


8.40(12) 


0.0664(14) 


So 


5.9138 


12 


4.601 


1.30 


10000 


2.501(36) 


0.0541(11) 


Bi 


6.0 


14 


5.368 


1.30 


10000 


2.761(40) 


0.0597(13) 


B 2 


6.0808 


16 


6.135 


1.30 


10000 


2.723(40) 


0.0589(13) 


Bz 


6.1568 


18 


6.902 


1.30 


10000 


2.788(40) 


0.0602(14) 



Table 1: Simulation parameters and results. The first seven lattices are gen- 
erated with the parameter value s = 0.4, while the other four, under the 
horizontal line, are for s = 0.0. The topological susceptibility is given by 
X = (Q^)/V = C 1 . 



centered around the point with L/a = 14. While we cannot exclude that 
the point at L/a = 12 turns out to be lower with respect to the others due 
to a statistical fluctuation only, the other three points indicate clearly that 
finite size effects are within our statistical errors for volumes larger or equal 
than (1.3 fm) 4 . 

Building on this result, the coupling constant of the other seven lattices 
have been chosen so that the lattice linear extent is fixed to be 1.3 fm, while 
the values of the coupling constant and of the negative mass s are chosen 
in order to properly estimate discretization errors at this level of precision. 
The results for the topological susceptibility for all lattices with linear ex- 
tent of 1.3 fm are shown in Fig. [2] as a function of (a/ro) 2 . The data show 
a clear trend to converge to the same value in the continuum limit within 
the statistical precision reached. Nevertheless a better understanding of the 
behaviour at small lattice spacing values, i.e. larger lattice volumes, seems 
to be necessary in order to improve the continuum limit study. In fact, 
data at the 2% level of error indicate a statistically non-negligible (and non- 
universal) contributions of terms of 0(a 4 ) to the behaviour of the two curves 
in this range of values of j3. Simulations at larger values of j3 are underway. 
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Figure 1: Rescaled topological charge as a function of the lattice size. Bands 
of ±2% and ±5% centered at the value measured at L = 14 are also shown. 
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Figure 2: Topological susceptibility as a function of the square of the lattice 
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